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Accurate accounting of particle number and 4-momentum in radiative transfer may be 
facilitated by the use of transport equations that allow transparent conversion between 
volume and surface integrals in both spacetime and momentum space. Such conservative 
formulations of general relativistic radiative transfer in multiple spatial dimensions are 
presented, and their relevance to core-collapse supernova simulations described. 

1. Neutrino transport in core-collapse supernovae 

Core-collapse supernovae result from the catastrophic collapse of the iron core of a 
massive star. Collapse is halted soon after the matter exceeds nuclear density. At 
this point a shock wave forms and begins to move out, but soon stalls as energy is 
lost to neutrino emission and dissociation of heavy nuclei falling through the shock. 
The details of how the stalled shock is revived sufficiently to continue plowing 
through the outer layers of the progenitor star are unclear. Some combination 
of neutrino heating of material behind the shock, convection, instability of the 
spherical accretion shock, rotation, and magnetic fields launches the explosion. 

It is natural to consider neutrino heating as a mechanism for shock revival. This 
is because neutrinos dominate the energetics of the post-bounce evolution, carrying 
away about 99% of the released gravitational energy. If we want to understand the 
explosion — which accounts for only about 1% of the energy budget of the system — 
we should carefully account for the neutrinos' much larger contribution to the energy 
budget, and accurately track both the energy and lepton number transferred from 
the neutrinos to the fluid. In principle, this requires tracking both the energy and 
angle dependence of the neutrino distribution functions at every point in space. A 
full treatment of this six-dimensional neutrino radiation hydrodynamics problem 
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remains too costly for currently available computational resources; the most that 
has been attempted so far has been a sort of "3.5" dimensional approximation^ 

2. The meaning of "conservative" 

Conservative formulations of radiation transport are one aspect of trying to get the 
details of neutrino radiative transfer right. To give an idea of the kind of formulation 
we seek, I first review "elemental" and "conservative" descriptions of the familiar 
dynamics of a fluid medium. 

The elemental formulation expresses the evolution of the fluid in terms of equa- 
tions of motion for the 4- velocity spatial components u 1 and two quantities measured 
by a comoving observer riding along with the fluid, the total energy density p and 
baryon number density n: 

-» dp ■ {P+P) 9 ' ^0=0, (2) 



dx^ \f—g dx^ 
J- (■f=gu fl )=0. (3) 



dx^ \J—g dx^ 

The name "elemental" denotes the fact that by writing down separate equations of 
motion for the velocity and two comoving-frame quantities, the kinetic and "intrin- 
sic" fluid energies — two "elements" of the system — are analytically separated. 

In contrast, the conservative approach expresses the evolution of the system in 
terms of the divergence of the stress-energy tensor T^ v = (p + p)u^u v + pg^ u and 
the divergence of a baryon number flux vector N 11 = nu^: 

' ° [V^gT^) = —T u pfl T pfi , (4) 



—g dx^ 

' d -(J—gN»)=Q. (5) 



\f—g dx^ 

Volume integrals of the left-hand sides — obtained by multiplying by the invariant 
spacetime volume element \/—g d x and integrating — are related to surface integrals 
in an obvious manner: 

1 d 



Physically, this formulation relates the time rates of change of 4-momentum and 
baryon number in a volume to fluxes through a surface surrounding that volume; 
hence the labeling of this formulation as "conservative." 

I will now describe why a common equation of particle (or radiation) transport — 
the Boltzmann equation — is an "elemental" formulation. The Boltzmann equation 
describes the evolution of /, the scalar distribution function: 

P^C^P- - P A ^V^^ - C[/]. (7) 
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The scalar distribution / is a function of time, spatial position, and momentum. 
Multiplied by suitable volume elements in phase space, it gives the number of par- 
ticles in these infinitesimal cells. The left-hand side of Eq. Q describes advection 
through phase space; it is closely related to the geodesic equations describing clas- 
sical particle trajectories. The right-hand side, the collision integral C, describes 
the change in the number of particles occupying a given trajectory. Note that the 
spacetime derivative is taken with respect to coordinates x M with a plain, unaccented 
index. This represents a global coordinate basis; in the spatially multidimensional 
case it will typically be a "lab frame" basis. A global time coordinate is convenient 
for numerically integrating the system forward in time. On the other hand, note the 
hats that accent the indices of the momenta. It is convenient to describe collisions 
in terms of momentum components measured with respect an orthonormal coordi- 
nate system carried by observers riding along with the fluid with which the neutrino 
radiation interacts. This "comoving basis" is indicated by the hatted indices. The 
tranformation £ M ^ relates the global coordinate basis to the comoving orthonormal 
basis, first through a transformation to a "lab frame" orthonormal basis, followed 
by a Lorentz boost to the comoving frame. Notice also that the momentum coor- 
dinates have been changed from p 1 to u 3 . This represents a change from Cartesian 
momentum coordinates to spherical momentum coordinates (energy and direction 
angles). Finally, note the connection coefficients T. These represent gravitational 
redshifts and trajectory bending induced by the curvature of spacetime, Doppler 
shifts and angular aberrations related to the Lorentz boost between the lab and 
comoving frames, and the use of curvilinear coordinate systems. 

In terms of the categories described above in connection with fluid evolution, 
the Boltzmann equation is "elemental" : The most fundamental quantity — the dis- 
tribution function — is the evolved variable, and volume integrals of the equation in 
both spacetime and momentum space are not obviously related to surface integrals. 

3. The conservative formulations 

Having discussed what I mean by "conservative," and having seen that the Boltz- 
mann equation is "elemental" or non-conservative, I now turn to conservative for- 
mulations of radiation transport,^ beginning with number conservation. The first 
momentum moment of the scalar distribution function constitutes a particle number 
flux vector: 



momentum bin to the lab frame number vector. The particle number flux vector 
N 11 satisfies a law of particle number conservation, 





where the quantity 



fiP^ f is the contribution of each comoving frame 
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which one would use to check number conservation. 

Noting that the right-hand side of Eq. © is simply the momentum integral 
of the collision integral, we see that this conservation law must simply be the mo- 
mentum integral of Eq. (JJJ). But it is not obvious how this works out, because of 
quantities outside the spacetime and momentum derivatives that depend on space- 
time position and momentum. Clearly, some conspiring cancellations are at work. 
But in the discretized case, in numerical work, a straightforward momentum inte- 
gration of a straightforwardly discretized Eq. Q will not be consistent numerically 
with a straightfoward discretization of Eq. @ . 

However, the following expression is equivalent to the Boltzmann equation, yet 
transparently related to Eq. @ : 

1 d 
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(10) 



Multiplying Eq. ^ by dP = e" 1 



du du du and integrating immedi- 



ately yields the number conservation law, Eq. @. Hence straightfoward differenc- 
ings of these two expressions will be numerically consistent. 

Another reformulation of the Boltzmann equation is an energy-conservative 
form, not detailed here. Suffice it to say that its momentum integral gives a trans- 
parent connection to the divergence of the particle stress-energy tensor. 

The number- and energy-conservative formulations — which respectively facil- 
itate computation of lepton number and energy transfer to essentially machine 
accuracy — might be used in a couple of different ways. Both formulations could 
be solved separately, in order to nail down the transfer of both energy and lepton 
number. The values of the scalar distribution function implied by these two differ- 
ent distributions would then serve as a consistency check. Alternatively, only one 
of the formulations might be solved, with the analytic relationship between the two 
conservative formulations being used to design finite difference expressions that 
provide consistency with the other formulation. This latter philosophy has been 
employed in the work of the group centered at Oak Ridge National Laboratory.^ 
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